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We applied multiresolution wavelet analysis to the sequence of times between human heartbeats 
(R-R intervals) and have found a scale window, between 16 and 32 heartbeats, over which the 
widths of the R-R wavelet coefficients fall into disjoint sets for normal and heart-failure patients. 
£ — , This has enabled us to correctly classify every patient in a standard data set as either belonging 

■ to the heart-failure or normal group with 100% accuracy, thereby providing a clinically significant 

^\ ' measure of the presence of heart-failure from the R-R intervals alone. Comparison is made with 

t-H , previous approaches, which have provided only statistically significant measures. 

PACS number(s): 87.10.+e, 87.80.+S, 87.90.+y 

►7- ■ Multiresolution wavelet analysis |j|-|5|] has proved to be a useful technique for analyzing signals at multiple scales, 
even in the presence of nonstationarities which often obscure such signals j^,|7j] . The sequence of times between human 
heartbeats (R-R intervals) is a prototype of a nonstationary time series that carries information about the state of 
' cardiovascular health of the patient [||| . 

By projecting this sequence into a wavelet space, a new set of variables is obtained, whose statistics allow us, for 
the first time, to correctly classify every patient in a standard data set as either heart-failure or normal, with 100% 
(yT) ', accuracy. It is clear from our results that the R-R intervals alone suffice as a measure for the presence of heart-failure; 
f^) ■ the full electrocardiogram in not required. This remarkable result arises from the ability of multiresolution analysis to 
simultaneously and compactly monitor multiple time scales and thereby to expose a hitherto unknown scale window 
(between 16 and 32 heartbeats) over which the widths of the R-R wavelet coefficients fall into disjoint sets for normal 
and heart-failure patients. The emergence of this particular scale window should help shed light on the underlying 
dynamics of cardiovascular function || . Previous approaches |ic] |lj] , even those that have made use of wavelets [[[3) , 
have been successful only in providing a statistically significant measure, rather than the clinically significant one we 
pOj have developed. The analysis method we have used is applicable to a wide variety of nonstationary physical and 
biological signals, regardless of whether the underlying fluctuations have stochastic origins or arise from nonlinear 
dynamical processes. 

The series of intervals between adjacent heartbeats t; (known as R-R or interbeat intervals in cardiology; see Fig. 
' la and lb) is thought to result from a complex superposition of multiple physiological processes at their respective 
characteristic time scales |J. The object of this letter is to demonstrate that is is possible, without any a priori 
knowledge of the physiological time scales or underlying heart dynamics, to determine a range of scales over which a 
. ! statistic of the wavelet coefficients permits each heart-failure and normal patient to be correctly categorized. 

Scale-dependent statistics are constructed by transforming the discrete-time sequence of R-R intervals s = {r^} 
fl4| into a space of wavelet coefficients. One can think of the transformed signal in terms of a landscape over a 
two-dimensional plane whose axes are interbeat-interval number i and scale m (see Fig lc). Smaller scales correspond 
to more rapid variations and therefore to higher frequencies. The height is the value of the corresponding wavelet 
coefficient. With such a three-dimensional construct it is possible to trace the relative importance of different scales as 
the heartbeat sequence proceeds. Technically the coefficients are obtained by carrying out a discrete wavelet transform 
(DWT) DH 

M 

^(«) = 2- m / 2 ^r^(2-^-n), (1) 
i=l 

where the scale variable m and the translation variable n are integers, and M represents the total number of R-R 
intervals analyzed. The discrete wavelet transform is evaluated at the points (m, n) in the scale-interval-number 
plane. 

We have carried out this transformation using a broad range of orthonormal, compactly supported analyzing 
wavelets. We present results for both Daubechies 10-tap and Haar wavelets; similar results were obtained using 
other wavelets. Orthogonality in the DWT provides that the information represented at a certain scale m is disjoint 
from the information at other scales. Because certain wavelets ip have vanishing moments, polynomial trends in the 
signal are automatically eliminated in the process of wavelet transformation |^,0,^5|. This is salutatory in the case 
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of the heartbeat time series, as is evident from the trends apparent in Fig. lb, which are eliminated by the wavelet 
transformation shown in Fig. Id. 

Since the signal s fluctuates in time, so too does the sequence of wavelet coefficients at any given scale, though 
its mean is zero [Q . In Fig. 2 we present histograms of the wavelet coefficients at six different scales (m = 2 - 7) 
for a normal (left) and a heart-failure patient (right). The wavelet coefficients for the heart-failure patient evidently 
exhibit substantially reduced variability, particularly at intermediate scales. A natural measure for this variability is 
the wavelet-coefficient standard deviation, as a function of scale: 



where N is the number of wavelet coefficients at a given scale m (N = M/2 m ). 

The principal results of this paper are displayed in Fig. 3, where cr wav is plotted vs. scale (1 < m < 10) for all of 
the 27 patients (open and solid circles) Oj, using the Haar wavelet (left) and the Daubechies 10-tap wavelet (right). 
At scales 4 and 5, corresponding to 2 4 -2° = 16-32 heartbeats, cr wav serves to completely separate the two classes 
of patients for both types of wavelets (white regions), thereby providing a clinically significant measure |^,^J of the 
presence of heart failure with 100% sensitivity at 100% specificity (such that all normals are so identified). One can 
do no better. Though it has been previously shown that complete separation can be achieved using heartrate analysis 
this is the first instance that we know of, in which the R-R intervals can be used as a definitive determinant 
of the presence of a heart disorder in an individual patient. Both at smaller, and at larger scales, there are multiple 
overlaps of the heart-failures and the normals, though the measure certainly remains statistically significant over all 
scales shown. Similar results are obtained for other analyzing wavelets. 

The results indicate that healthy patients exhibit greater fluctuations than those afflicted with heart failure over 
a time scale of 16-32 heartbeats (roughly 0.2 to 0.5 minutes). This also appears to apply for sudden cardiac death 
(SCD) as illustrated by the open squares in Fig. 3. It is tempting to ascribe the physiological origin of this window 
to baroreflex modulations of the sympathetic or parasympathetic tone, which lie in the range 0.04 to 0.09 Hz (0.2 to 
0.5 minutes), but we do not believe that this is correct. Rather, we expect that this window likely has its origin in the 
intrinsic behavior of the heart itself. It will be important to carry out a thorough study in which our multiresolution 
wavelet-analysis technique is applied to the R-R intervals from transplanted hearts to assess the role that the autonomic 
system might play in heartrate variability. 

It is useful to tease apart the roles played by the magnitudes Ti of the interbeat intervals and their ordering in 
achieving this complete separation. The effects of the former continue to reside in the randomly reordered (shuffled) 
sequence of R-R intervals; however information about the ordering is removed in this surrogate data set p7| . We 
therefore calculate the standard deviation <J^" £ for all 27 heartbeat time series after shuffling the R-R intervals. 
[The first 70 000 interbeat intervals were selected for analysis after the entire data sets (see Table 1 in Ref. [pT[) 
were shuffled.] The results are shown in Fig. 4a. It is clear that the two classes of patients are no longer completely 
separated; three heart-failure patients fall among the normals at all scales, yielding a sensitivity of 80% at a specificity 
of 100%. Comparison with Fig. 4b shows that the shuffled- wavelet result is essentially identical to that obtained by 
using the standard deviation cr; n t of the interbeat-interval histogram (IIH) obtained from these data sets, a measure 
that has long been used in cardiology p9| , p0[ . The concurrence in sensitivity displayed in Figs. 4a and 4b is not 
accidental; the shuffled interbeat intervals essentially comprise a renewal process so that the IIH contains all of the 
available information. All dependencies among intervals, and therefore long-term correlations, are removed from the 
shuffled surrogate data jL7|, leaving behind only short-term information. Indeed, for the Haar analyzing wavelet, 
a wa¥ (m = 0) (in the absence or in the presence of shuffling) is analytically identical to <7i n t. 

The ordering of the interbeat intervals gives rise to scaling behavior, as is evident from a comparison of Figs. 3 
and 4a. For normal patients (open circles, solid lines) the straight-line behavior in Fig. 3 indicates that approximate 
scaling is maintained across all scales whereas for heart-failure patients (filled circles, dashed lines) the relatively flat 
nature of the curves in the region m < 3 indicates that cr wav is essentially scale independent in this region. 

The distinction can be examined quantitatively by calculating the average scaling exponents a |2l|] in the two 
ranges (1 < m < 3 and 3 < m < 10), for both classes of data. The results are reported in Table 1, along with the 
average scaling exponents determined across the entire scaling region examined (1 < m < 10). It is clear that for the 
12 normal patients, the value of a is quite insensitive to the range over which it is estimated; it is nearly the same 
in the two subregions as it is in the entire range. For the 15 heart-failure patients, on the other hand, the scaling 
exponent estimated in the region 1 < m < 3 is dramatically lower than that estimated in the region 3 < m < 10. 
Over the range of larger scales, the values of the scaling exponents a provided by our wavelet measure (see Table 1), 
are in good accord with typical values S provided by the interval-based periodogram at low frequencies (see Table 1 
in Ref. 0|), as expected. We conclude that scaling persists across a broader range for normal patients than it does 
for heart-failures, as is visually evident in Fig. 3. 
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These observations lead us to consider a heart-failure index determined by the difference of these scaling exponents: 
A = a(3 < m < 10) — a(l < m < 3). Evaluating A we find that two heart-failures fall among the normals, 
corresponding to a sensitivity of 87% at a specificity of 100%. Thus considering only the scaling information in cr wav , 
while ignoring the magnitude differences for normals and heart-failures associated with short-term information (as 
illustrated in Fig. 4a), fails to give rise to complete separation. 

Over the years, using this same collection of data, a number of measures based on scaling have been evaluated for 
their accuracy in discriminating between normal and heart-failure patients. Peng et al. JToj] examined the correlation 
properties of the heartbeat-interval increments Ii — {t^+i — Tj}, and obtained the exponent of the associated power-law 
spectrum, which they denoted as @. It was shown subsequently JL7| that this measure is isomorphic to the exponent 5 
of the interval-based spectrum [£2| at low frequencies |2^] and therefore reveals only long-term correlations. We have 
calculated f3 for all 27 data sets and present the results in Fig. 4c. This measure results in 7 heart-failures among the 
normals, yielding a sensitivity of 53% at 100% specificity, so that it is not particularly successful in discriminating the 
two classes of patients. The interbeat-interval standard deviation cr int shown in Fig. 4b (sensitivity of 80% at 100% 
specificity) performs significantly better. In a subsequent paper, Peng and collaborators [jllj constructed a "detrended 
fluctuation measure" which turns out to yield approximately 80% sensitivity at 100% specificity for the 27 data sets. 
Most recently, this same group introduced a so-called "scaling instability index" fL2|| , which achieved a sensitivity 
of 71% at 100% specificity for the 25 (of the same 27) data sets that they reported, as shown in Fig. 4d. The 
performance of both of these latter measures is therefore comparable to that achievable with the interbeat-interval 
standard deviation measure cr int introduced by Wolf et al. long ago [fl9|| . As far as we are aware, no group has been 
able to successfully separate this standard collection of data with 100% sensitivity at 100% specificity using scaling 
measures alone. 

We conclude that our multiresolution approach succeeds not only because it eliminates trends in a mathematically 
acceptable way, but also because it crisply reveals a range of scales over which heart-failure patients differ from 
normals, both in short- and long-term heartbeat behavior. In contrast, interbeat-interval measures reflect only short- 
term behavior, whereas scaling measures reflect only long-term behavior. Our approach should have wide applicability 
in the analysis of nonstationary processes in the physical and biological sciences. 

We are grateful to S.B. Lowen for valuable comments. 



FIG. 1. (a) Schematic diagram of an electrocardiogram segment, showing the beat occurrence times ti and the interbeat 
(R-R) intervals t%. (b) Series of interbeat intervals Ti versus interval number i for a typical normal patient (data set 16265). 
(Adjacent values of the interbeat interval are connected by straight lines to facilitate viewing.) Substantial trends are evident, 
(c) 3-D representation of the wavelet coefficient W as a function of scale (1 < m < 10) and interval number i, over a portion 
of the data set, using a Daubechies 10-tap analyzing wavelet, (d) Wavelet coefficient at three scales (m = 2, 4, 8) for the data 
set illustrated in Fig lb. The trends in the original interbeat-interval time series are removed by the wavelet transformation. 



FIG. 2. Typical histograms of the interbeat-interval wavelet coefficients for a normal patient (data set 16265, left panels) 
and for a heart-failure patient (data set 6796, right panels), for different values of the scale m, using a Daubechies 10-tap 
wavelet. The ordinate is the number of wavelet coefficients N. Heart-failure patients tend to have narrower distributions, 
particularly at intermediate scales. 



FIG. 3. Wavelet-coefficient standard deviation <r wav versus scale m for the standard 27 data-set collection (12 normals and 
15 heart-failures indicated by circles) [14], using the Haar wavelet (left) and the Daubechies 10-tap wavelet (right). Complete 
separation of the two groups is achieved at scales 4 and 5, corresponding to 2 4 - 2 5 heartbeats. Results for an SCD patient 
(white squares), using the same number of interbeat intervals, also exhibit low variability. The outcome is similar for both 
analyzing wavelets, though the separation is slightly more pronounced for the Daubechies wavelet, most likely because of its 
greater number of vanishing moments. 
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FIG. 4. (a) Wavelet-coefficient standard deviation a^, 1 versus scale m for the standard 27 data-set collection [14] after 
random reordering of the R-R intervals (shuffling), using a Daubechies 10-tap wavelet. Partial separation of the two groups is 
achieved; three heart-failure patients fall among the normals at all scales, yielding a sensitivity of 80% at a specificity of 100%. 
(b) Standard deviation of the interbeat intervals a- lnt for all 27 data sets. This measure was first used by Wolf et al. [19]. 
Three heart-failure patients fall among the normals, corresponding to a sensitivity of 80% at 100% specificity, (c) Exponent 
f3 of the heartbeat-interval increment-process spectrum for all 27 data sets; this measure was used by Peng et al. [10]. Seven 
heart-failures fall among the normals, yielding 53% sensitivity at 100% specificity, (d) The "scaling instability index," recently 
introduced by Peng and collaborators [12], and reported by them for 25 of the same 27 data sets. Four heart-failures fall among 
the normals, yielding 71% sensitivity at 100% specificity, still falling short of our multiresolution results shown in Fig. 3 which 
yield 100% sensitivity at 100% specificity. 
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a(l < m < 10) 


a(l < m < 3) 


a(3 < m < 10) 


Normal (12) 


0.79 ± 0.08 


1.23 ±0.13 


1.40 ±0.37 


1.22 ±0.11 


HF (15) 


0.67 ±0.13 


1.35 ±0.22 


0.26 ±0.60 


1.57 ±0.17 



TABLE I. Mean interbeat intervals and estimated scaling exponents a (with their standard deviations) for the 12 normal 
and 15 heart-failure patients, collected in two separate groups. Scaling-exponent estimates are evaluated over the whole range of 
m, as well as over two smaller scaling subregions: 1 < m < 3 and 3 < m < 10. The scaling behavior is significantly different for 
the two groups of patients. The difference of the scaling exponents in the two subregions, A = a(3 < m < 10) — q(1 < m < 3), 
leads to a sensitivity of 87% at 100% specificity. 
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